arXiv:1501.00252vl [astro-ph.SR] 1 Jan 2015 


Astronomy & Astrophysics manuscript no. MSMSK'rev'finar4a @ ESO 2015 

January 5, 2015 


Torsional Alfven Waves in Solar Magnetic Flux Tubes of Axial 

Symmetry 

K. Murawski^, A. Solov’ev^, Z.E. Musielak^’^, A.K. Srivastava^, and J. Kraskiewicz^ 

^ Group of Astrophysics, University of Maria Curie-Sklodowska, ul. Radziszewskiego 10, 20-031 Lublin, Poland 
e-mail: kmur@kft.umcs.lublin.pl 

^ Central (Pulkovo) Astronomical Observatory, Russian Academy of Sciences, St. Petersburg, Russia 
e-mail: solov@gao.spb.ru 

^ Department of Physics, University of Texas at Arlington, Arlington, TX 76019, USA 
e-mail: zmusielak@uta.edu 

^ Kiepenheuer-Institut fiir Sonnenphysik, Schoneckstr. 6, 79104 Freiburg, Germany 
e-mail: musielak@kis . uni-f reiburg. de 

^ Department of Physics, Indian Institute of Technology (Banaras Hindu University), Varanasi-221005, India 
e-mail: asrivastava. app@iitbhu. ac. in 

Received jdate^ / Accepted jdate^ 


ABSTRACT 

Aims. Propagation and energy transfer of torsional Alfven waves in solar magnetic flux tubes of axial symmetry is 
studied. 

Methods. An analytical model of a solar magnetic flux tube of axial symmetry is developed by specifying a magnetic 
flux and deriving general analytical formulae for the equilibrium mass density and a gas pressure. The main advantage 
of this model is that it can be easily adopted to any axisymmetric magnetic structure. The model is used to simulate 
numerically the propagation of nonlinear Alfven waves in such 2D flux tubes of axial symmetry embedded in the solar 
atmosphere. The waves are excited by a localized pulse in the azimuthal component of velocity and launched at the 
top of the solar photosphere, and they propagate through the solar chromosphere, transition region, and into the solar 
corona. 

Results. The results of our numerical simulations reveal a complex scenario of twisted magnetic field lines and flows 
associated with torsional Alfven waves as well as energy transfer to the magnetoacoustic waves that are triggered by 
the Alfven waves and are akin to the vertical jet flows. Alfven waves experience about 5% amplitude reflection at the 
transition region. Magnetic (velocity) field perturbations experience attenuation (growth) with height is agreement with 
analytical findings. Kinetic energy of magnetoacoustic waves consists of 25% of the total energy of Alfven waves. The 
energy transfer may lead to localized mass transport in the form of vertical jets, as well as to localized heating as slow 
magnetoacoustic waves are prone to dissipation in the inner corona. 
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1. Introduction 

The role of magnetohydro dynamic (MHD) waves in trans¬ 
porting energy and heating various layers of the solar atmo¬ 
sphere, and acceleration of the solar wind, has been studied 
by many authors (e.g., Hollweg 1981; Hollweg et al. 1982; 
Priest 1982; Hollweg 1992; Musielak 1998; Ulmschneider 
& Musielak 1998, 2003; Dwivedi & Srivastava 2006, 2010; 
Zaqarashvili & Murawski 2007; Gruszecki et al. 2008; 
Murawski & Musielak 2010; Ofman 2009; Chmielewski et 
al. 2013). Among the three basic MHD modes, Alfven waves 
are of a particular interest as they can carry their energy 
along magnetic field lines to higher atmospheric altitudes 
on fast temporal scales (e.g., Hollweg 1978, 1981, 1992; 
Musielak & Moore 1995; Cargill et al. 1997; Vasheghani 
Farahani et al. 2011, 2012). 

The interest in Alfven waves has recently significantly 
increased because of their observational evidence in the so¬ 
lar atmosphere. Several authors claim that they have al¬ 


ready observed outwardly-propagating Alfven waves in the 
quiescent solar atmosphere in the localized solar magnetic 
fiux tubes (e.g., Jess et al. 2009; Fujimura & Tsuneta 2009; 
Okamoto & De Pontieu 2011; De Pontieu et al. 2012; Sekse 
et al. 2013), while Van Doorsselaere et al. (2008) questioned 
the interpretation of Alfven wave observations in the solar 
atmosphere. Specifically, Jess et al. (2009) interpreted the 
Ha data in terms of torsional Alfven waves in the solar chro¬ 
mosphere, with periods ranging from 2 min to near 12 min, 
and with the maximum power near 6 — 7 min. Torsional 
and swirl-like motions were also reported by respectively 
Bonet et al. (2008) and Wedemeyer-Bohm & Rouppe van 
der Voort (2009). According to Bonet et al. (2008), there are 
vortex motions of G band bright points around downflow 
zones in the photosphere, and lifetimes of these motions 
are of the order of 5 min. Wedemeyer-Bohm & Rouppe van 
der Voort (2009) observed disorganized relative motions of 
photospheric bright points and concluded that they induce 
swirl-like motions in the solar chromosphere. 
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Different theoretical aspects of the generation, propa¬ 
gation and dissipation of torsional Alfven waves in the so¬ 
lar atmosphere were investigated by Hollweg (1978, 1981), 
Heinemann & Olbert (1980), and Ferriz-Mas et al. (1989), 
who solved the Alfven wave equations using different sets 
of wave variables. Musielak et al. (2007) demonstrated that 
the propagation of these waves along isothermal and thin 
magnetic flux tubes is cutoff-free. However, Routh et al. 
(2007, 2010) showed that gradients of physical parameters, 
such as temperature, are responsible for the origin of cutoff 
frequencies, which restrict the wave propagation to certain 
frequency intervals. In most of the above work only linear 
(small-amplitude) Alfven waves were considered. 

Numerical studies of linear and nonlinear torsional 
Alfven waves in solar magnetic flux tubes were performed 
by Hollweg et al. (1982), who considered a one-dimensional 
(ID) model. The original work of Holweg et al. was ex¬ 
tended to higher dimensions by Kudoh & Shibata (1999) 
and Saito et al. (2001), who investigated the effects caused 
by nonlinearities. Fedun et al. (2011) showed by means of 
numerical simulations that chromospheric magnetic flux- 
tubes can act as a frequency filter for torsional Alfven 
waves. Wedemeyer-Bohm et al. (2012) used the observa¬ 
tional results originally reported by Wedemeyer-Bohm & 
Rouppe van der Voort (2009) to demonstrate that the swirl¬ 
like motions observed in the lower parts of the solar at¬ 
mosphere produce magnetic tornado-like motions in the 
solar transition region and corona. However, more recent 
studies performed by Shelyag et al. (2013) seem to imply 
that the tornado-like motions actually do not exist once 
time dependence of the local velocity field is taken into 
account. Instead, the tornado-like motions were identified 
as torsional Alfven waves propagating along solar mag¬ 
netic flux tubes. Recently, Wedemeyer-Bohm & Rouppe 
van der Voort (2009) reported on small, rotating swirls ob¬ 
served in the solar chromosphere, and interpreted them as 
plasma spiraling upwards in a funnel-like magnetic struc¬ 
ture. Chmielewski et al. (2014) and Murawski et al. (2014) 
modeled such short-lived swirls and associated plasma and 
torsional motions in coronal magnetic structures. 

As a result of the expanding, curved and strong mag¬ 
netic field lines of solar magnetic flux tubes embedded in the 
solar atmosphere, formulating a model of such flux tubes is 
always a formidable task. Nevertheless, many efforts were 
undertaken in the past to develop realistic flux tube models. 
Among many others, let us mention Low (1980) and Gent 
et al. (2013, and references therein), who constructed mag¬ 
netic flux tube models in the solar atmosphere. In partic¬ 
ular, Low (1980) formulated an analytical sunspot model, 
while Gent et al. (2013) solved analytically the magnetohy¬ 
drostatic equilibrium problem. In this paper, we develop a 
model of a solar magnetic flux tube of axial symmetry. The 
model requires specifying a magnetic flux and gives ana¬ 
lytical formulae for the equilibrium mass density and a gas 
pressure; these formulae are general enough to be applied 
to any axisymmetric magnetic structure. 

The main goal of this paper is to perform numerical 
simulations of nonlinear torsional Alfven waves in our de¬ 
veloped model of magnetic flux of axial symmetry tube 
embedded in the solar atmosphere with the temperature 
profile given by Avrett & Loeser (2008). We also study the 
energy transfer from Alfven waves to slow magnetoacous¬ 
tic waves in the non-linear regime, as well as their reflec¬ 
tion and transmission through realistic solar transition re- 



y [7l/m] 

Fig. 1. Hydrostatic equilibrium profile of solar temperature. 



Fig. 2. Equilibrium magnetic field lines. Red, green, and blue 
arrows correspond to the x-, y-, and z-axis, respectively. The size 
of the box shown is (—0.75,0.75) x (0,4) x (—0.75,0.75) Mm. 
The color map corresponds to the magnitude of a magnetic field. 

gion, which are significant in constituting vertical jet flows 
(mass) and localized energy transport processes in the in¬ 
ner corona. The paper is organized as follows. Our model of 
a solar magnetic flux tube of axial symmetry embedded in 
the solar atmosphere is developed in Section 2. The results 
of our numerical simulations are presented and discussed in 
Section 3. Gonclusions are given in Section 4. 

2. A model of a solar magnetic flux tube of axial 
symmetry 

We consider a model solar atmosphere that is described by 
the following ideal, adiabatic, 3D magnetohydro dynamic 
(MHD) equations: 

|f + v-(W) = o, (1) 
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Fig. 3. Equilibrium profiles of the Alfven (solid) and sound (dashed) speeds (the left panel) and plasma jS (the right panel) along 
the flux tube axis. 



X (Mm) 


4.8 



1-9.6 

(km s '!) 


Fig. 4. Numerical blocks with their boundaries (solid lines) and 
the pulse in velocity Vz of Eq. (l3Ql) (color maps) in the x — y 
plane for z = 0 Mm at t = 0 seconds. A part of the simulation 
region is only displayed. 


d\ 


-( 


dp 

dt 


dB 

dt 


') X B + £»g , 

(2) 

p = — qT, 

(3) 

m 

V-B = 0 , 

(4) 


We assume that the solar atmosphere is in static equi¬ 
librium (Ve = 0) with the Lorentz force balanced by the 
gravity force and the gas pressure gradient, which means 
that 

-(V X Be) X Be + Qeg “ Vpe = 0 , (5) 

P 

where the subscript 'e' corresponds to the equilibrium con¬ 
figuration. We consider an axisymmetric flux tube whose 
equilibrium is described by Eq. (|5j) and its magnetic field 
satisfies the solenoidal condition. We now introduce the 
magnetic fiux (^) normalized by 27r, and obtain 

[ B^yv'dr', ( 6 ) 

Jo 

where Bey is a vertical magnetic field component and r is 
the radial distance from the flux tube axis given by 


r = + 2 )^ . (7) 

As a result of Eq. the magnetic field is automatically 
divergence-free and its radial (5er), azimuthal (Bee) and 
vertical (Bey) components are expressed as 


Be 


R -n B -1^ 

rA 5 ^eO — b , iDey — oj 

r dy r dr 


(8) 


We consider a magnetic fiux tube of axial symmetry, which 
is initially non-twisted, and represent the tube magnetic 
field in terms of its azimuthal component of the vector- 
potential (AO) as 


Be = V X A<9 , 


(9) 


where ^ is a unit vector along the azimuthal direction. In 
this case, we have 


where g is mass density, V is the flow velocity, and B 
is the magnetic field. The standard notation is used for 
the other physical parameters, namely, p, T, 7 = 5/3, /r, 
g = ( 0 , —^, 0 ), m, and /cp are the gas pressure, tempera¬ 
ture, adiabatic index, magnetic permeability, mean particle 
mass, and Boltzmann’s constant, respectively. The magni¬ 
tude of the gravitational acceleration is ^ = 274 m s“^. 
The value of m is specified by the mean molecular weight, 
which is 1.24 in the solar photosphere (Oskar Steiner, pri¬ 
vate communication) and assumed to be constant in the 
entire flux tube model. 


dA 

Bqv — 5 BeO — 0 , 

dy 


_ 1 d{rA) 


( 10 ) 


Comparing Eqs. (0 and m, we find 


Mr,y) = -'^{r,y). (11) 

r 

In Cartesian coordinates the magnetic vector potential 
components are 


[Ax, Ay, Az] 


r . ^ ^ .X, 

A-, 0, —A— . 

r r 


( 12 ) 
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These formulae are useful for uumerical simulatious because 
they ideutically satisfy the seleuoidal couditiou, which is 
importaut iu coutrolliug uumerical errors. 

Multiplyiug uow Eq. by Be, we obtaiu 


Be • (Vpe - Qeg) = 0 . (13) 

Since Pe = u), therefore, from the above equation with 

the use of Eq. m, we have 


dPe{^,y) 

=-b- 

dy 


(14) 


which meaus that the hydrostatic couditiou is satisfied 
aloug maguetic field hues, aud that ^ = coust. 

Rewritiug Eq. (|5|) iu its compoueuts aud takiug Eq. (|8|) 
iuto accouut, we obtaiu the so-called Grad-Shafrauov equa- 
tiou (e.g.. Low 1975; Priest 1982) giveu by 


However, it follows from Eq. ([Tb]) that we ueed to specify 
S{r,y) as 


- 2 r^\dr) ^ J ^ 2/2 dr r2 ' 


( 21 ) 


We calculate first dS{r,y)/dr, and then use Eq. to 
determine dS{y, ^)/d^. The result is 


dpe{y,'^) 


1_5 
r dr 



1 g 2 ^- 

y,2 Qy2 


( 22 ) 


Usiug Eq. we get 


dpejy,'^) ^ dpe{r, y) _ dp^{y,'^) d^ 

dy dy d^ dy 

Eiually, with help of Eq. (|2^ . we write 


^ ^ = ,:r.2 dPe{^,y) 

dr2 r dr dy'^ ^ d^ ' ^ ^ 

If we cousider the equilibrium gas pressure as a giveu 
fuuctiou, theu Eq. m becomes a uouliuear Dirichlet prob¬ 
lem for Iu coutrast. Low (1980) proposed the iuverse 
problem by specifyiug maguetic field through a choice of ^ 
by iutegratiug Eq. aud deriviug the geueral formulas 
for Pq aud Qq giveu iu terms of ^ aud its derivatives. Iu our 
approach, we regard ^ as a parameter, aud after iutegratiug 
Eq. (HU over we hud (Solov’ev 2010) 


yy dy 


1 



with 

Qh{y) 


Phjy) 

gMy) 


0 e{r,y) = Qh{y) + 
fd^'^d^dr 

\dy ) j ^ J dy"^ dr r"^ 


Id^ d 

r dy dr 



},(24) 


(25) 


Peir,y) =Phiy) - 

/i 


1 /d^y fdHd^dr 
2 r2 \ dr J J dy"^ dr r2 


, (16) 


where 


Ph{y) =Po exp 


r dy' 

Jy, Hy') 


(17) 


beiug the hydrostatic mass deusity. 

Note that iu the above formulas the maguetic fiux (^) 
aud cousequeutly the maguetic flux fuuctiou {A) are free to 
choose. Eor a flux tube, we cau specify them as 

A{r,y) = Boexp(-fc22/2)__L_ + ^Byor , (26) 

^{r,y) =rA{r,y), (27) 


is the hydrostatic gas pressure, aud 


A(y) = MkM 

mg 


(18) 


is the pressure scale-height aud T\^{y) is a hydrostatic tem¬ 
perature profile. 

We adopt T\^{y) for the solar atmosphere that is speci¬ 
fied by the model developed by Avrett & Loeser (2008). 
This temperature profile is smoothly exteuded iuto the 
coroua (Eig. [T]). It should be uoted that iu our model the 
solar photosphere occupies the regiou 0 < ^ < 0.5 Mm, the 
chromosphere is saudwiched betweeu ^ = 0.5 Mm aud the 
trausitiou regiou, which is located at ^ ^ 2.1 Mm, aud that 
above this height the atmospheric layers represeut the solar 
coroua. 

Accordiug to Eq. m, the equilibrium mass deusity 
cau be calculated if dpe{y, '^)/dy is kuowu. Moreover, from 
Eq. (p!6|) . we get Pe = p^ipr^y). Thus, iu order to hud 
dpe{y, ^)/5^, we use the followiug relatious that are valid 
for auy differeutiable fuuctiou (S): 


dS{r,y) 

dr 

dS{r,y) 

dy 


dS{^,y) 
d^f dr 

dS{^,y) dS{^,y) d^ 

dy d^ dy 


(19) 

( 20 ) 


where Hyo is the exterual maguetic field aloug the vertical 
directiou, aud ky are iuverse leugth scales aloug the ra¬ 
dial aud vertical directious, respectively. The maguetic field 
hues, which follow from Eqs. m aud (j26j) are displayed iu 
Eig. [2] for kr = ky = 4 Mm“^. We choose the maguitude of 
the refereuce maguetic field Bq iu such way that the mag¬ 
uetic field withiu the flux tube, at {x = O^y = 0.5,2) = 0) 
Mm, is about 285 Gauss, aud Byo ^ 11A Gauss. Eor these 
values, the resultiug maguetic field hues are predomiuautly 
vertical arouud the hue, x = z = 0 Mm, while further out 
they are beut aud Be decays with distauce from this hue. 
It must be also uoted that the maguetic field at the top of 
the simulatiou regiou is esseutially uuiform, with its value 
of Byp. 

Eigure [3] illustrates the souud speed, Cg, (the dashed hue 
iu the left pauel), Alfveu speed, ca, (the solid hue iu the 
same pauel), aud plasma /3 (iu the right pauel) giveu by 


Csir,y) = 


hpejr, y) 

Qe{r,y) 


, CA{r,y) = 


P{r, y) = 


/ Bl{r,y) 
y0e{r,y) '' 
2 ype{r, y) 
Bl{r, y) 


(28) 

(29) 


The souud speed below the trausitiou regiou is smaller thau 
10 km s“^, however, it abruptly raises iu the solar trausitiou 




































K. Murawski et aL: Non-linear Alfven Waves in Magnetic Flux Tubes 


5 


region, and then reaches its coronal value of 100 km s“^ at 
the height ^ = 10 Mm (not shown in the figure). The Alfven 
speed reveals similar trends, reaching a value of more than 
700 km s“^ at ^ = 10 Mm (also not shown). The coronal 
plasma /3is~3x 10“^ at ^ = 4 Mm (the right panel). It 
increases with the depth but remains lower than 1 above 
the height ^ ^ 0.8 Mm. It should be noted that as a result 
of strong magnetic field, /3 is lower within the fiux tube 
than in the surrounding medium (not shown). 


3. Results of numerical simulations 


Equations 0-0 are solved numerically using the code 
FLASH (Lee & Deane 2009; Lee 2013). This code imple¬ 
ments a third-order unsplit Godunov solver (Murawski & 
Lee 2011) with various slope limiters and Riemann solvers, 
as well as adaptive mesh refinement. We choose the van Leer 
slope limiter and the Roe Riemann solver (Murawski & Lee 
2012). For all the considered cases, we set the simulation 
box as (—1.25,1.25) Mm x (0, 20) x (—1.25,1.25) Mm, and 
impose boundary conditions by fixing in time all plasma 
quantities at all six boundary surfaces to their equilibrium 
values. In the present numerical studies, we use a non- 
uniform grid with a minimum (maximum) level of refine¬ 
ment set to 2 (5). The grid system at t = 0 s is shown in 
Fig.lUwith the blocks being displayed only up to ^ = 6 Mm. 
Above this altitude the block system remains homogeneous. 
As each block consists of 8 x 8 x 8 identical numerical cells, 
we reach the effective finest spatial resolution of 19.53 km, 
below the altitude ^ = 4.25 Mm. 

We initially perturb the azimuthal component of veloc¬ 
ity, Vo, by using 


V0e=[V^,V,] = ^vexp 

w w 


r'^ + {y - yof 

w‘^ 


, (30) 


where denotes the amplitude of the pulse, 0 is the unit 
vector along the azimuthal direction, yo is its vertical posi¬ 
tion, and w is the width. We set = 150 km s“^, w = 150 
km, and yo = 500 km, and hold them fixed. This value of 
Ay results in the effective maximum velocity of about 9.6 
km s“^ (Fig. 0]). 

The initial velocity pulse of Eq. OT triggers torsional 
Alfven waves. These waves are illustrated in Fig. O which 
shows magnetic field lines at t = 50 s (the top-left panel), 
t = 150 s (the top-right panel), t = 200 s (the bottom-left 
panel), and t = 250 s (the bottom-right panel). At t = 
50 s (the top-left panel) the upwardly propagating Alfven 
waves, while observed in Bq, reach the altitude y ^ I Mm, 
and at a later time Alfven waves spread in space while 
propagating along diverged magnetic field lines (the top- 
right and bottom panels). 

Note that small-amplitude magnetic field and velocity 
perturbations that include the effect of magnetic field ex¬ 
pansion, are governed by the following wave equations (Eqs. 
(3) & (4) in Hollweg 1981): 


(r^f) _ r 2^ 

dt^ j-lQe ds ^ ® ’ 

d% , d ( B, dh\ 
dt^ ^ds\^iQ^r^ ds) ' 


(31) 

(32) 


where v = Vq/ r, b = vBq, Vq and Bq are the azimuthal 
components of respectively perturbed velocity and mag¬ 


netic field, and s is the coordinate along the magnetic field 
lines. 

Moran (2001) showed that Ve oc Qe^^^jAB^, where A 
is the cross-sectional area of the flux tube. Near the flux 
tube axis we can take ABq approximately constant, hence 
in that region Vq (x On the other hand, Bq cx 

(Moran 2001). Hence, near the flux tube axis, these quan¬ 
tities are only dependent on the variation of Qq, with Vq 
and Bq respectively growing and decreasing with height. 
As a result of that the perturbations of magnetic field lines 
remain weak in the solar corona (see also Murawski & 
Musielak 2010 for the corresponding analysis in the case 
of the straight vertical magnetic field lines). 

Ratio of magnetic to kinetic energies (evaluated in the 
entire computational domain) of Alfven waves is shown in 
Fig. [7] (top). Initially, at t = 0 s, as we launch the initial 
pulse in the azimuthal velocity alone the magnetic energy 
is zero at t = 0 s but, according to our numerical simula¬ 
tions, contributions of the magnetic energy grow in time. It 
is well known that for torsional Alfven waves, there is the 
equipartition between potential and kinetic energies. Hence 
according to the top panel of Fig. [71 this equipartition of en¬ 
ergy is not reached in the numerical simulation until about 
t = 200 s (where the energy ratio equals unity). This is the 
natural explanation for the limiting value of the energy ra¬ 
tio being approximately unity at t = 200 s. At t = 225 s the 
magnetic energy is higher than the kinetic energy with the 
ratio equal to ~ 1.1. Ratios of slow magnetoacoustic waves 
energy to Alfven waves energy (asterisks) and fast magne¬ 
toacoustic waves energy to Alfven waves energy (diamonds) 
is illustrated in Fig. [7] (bottom). At t = 0 s kinetic energies 
of both the slow and fast magnetoacoustic waves are zero 
but as these waves are driven by the ponderomotive force 
which results from Alfven waves the energy ratio grows in 
time reaching a maximum level of about 0.25 which means 
that about 25% of the Alfven wave energy was converted 
into kinetic energies of magnetoacoustic waves. These waves 
propagate along magnetic field lines and power the vertical 
plasma flows in form of the jets. 

Figure [8] illustrates maximum of Vq (left-panel) and Bq 
(right-panel) as evaluated from the numerical data. We 
clearly see that Vq grows with height as predicted on the 
basis of the liner theory approximation reported by Moran 
(2001). On the other hand, Bq exhibits a fall-off with y with 
the numerical data points being close to the linear theory 
data (solid-line) as reported by Moran (2001). It must also 
be noted that the relationships shown by Moran (2001) are 
only valid when there is the equipartition of energy for tor¬ 
sional Alfven waves. According to the top panel of Fig. [71 
this occurs only the time of our numerical simulations is 
t = 200 s or longer. As the results of Fig. [ 8 ] are given for 
the times shorter than t = 200 s, then comparison of these 
results with the analytical findings of Moran (2001) is not 
strictly valid. This fall-off is also clearly evident in the spa¬ 
tial profiles of Bz{x, y,z = 0) as displayed at t = 150 s and 
t = 200 s (Fig. [ 6 l bottom panels). At t = 50 s and t = 100 
s (top panels), the perturbations of B^{x,y,z = 0) are lo¬ 
calized below the transition region but at later times the 
perturbations that penetrate the inner solar corona are of 
significantly lower magnitude than those remaining below 
the transition region. Consequently, these profiles demon¬ 
strate clearly that Alfven waves, while traced in perturbed 
magnetic field lines, are essentially located below the transi- 
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Fig. 5. Magnetic field lines at t = 50 s (top-left), t = 150 s (top-right), t = 200 s (bottom-left), and t = 250 s (bottom- 
right). Red, green and blue arrows design the directions of the x-, y-, and z-axis, respectively. The size of the box shown is 
(—0.75,0.75) X (0,4) x (—0.75,0.75) Mm. The color map corresponds to the magnitude of a magnetic held. 


tion region, which agrees with the conclusions of Murawski 
& Musielak (2010) who used the analytical and numerical 
techniques for the ID wave equations and found that mag¬ 
netic held perturbations, while excited below the transition 
region, remain week in the solar corona in comparison to 
the perturbations in lower atmospheric layers. 

Figure [9] shows temporal snapshots of velocity stream¬ 
lines at t = 100 s (the left-top panel), t = 150 s (the right- 
top panel), t = 200 s (the left-bottom panel), and t = 225 s 
(the right-bottom panel). These streamlines are dehned by 
the following equations: 

dx dy dz , . 

Note that a signal in the azimuthal velocity is clearly seen in 
the solar corona. Indeed, the streamlines show that velocity 


perturbations penetrate the transition region and enter the 
solar corona, leading to a generation of the main centrally 
located vortex, which starts developing at t = 100 s. 

Our results clearly show that the velocity perturba¬ 
tions are significantly enhanced in the solar transition re¬ 
gion as a result of very steep temperature gradient. The 
main physical consequence of this decrease is that in the 
inner corona, just above the transition region, we see the 
more velocity streamline perturbations (Figs. [9]), while we 
do not see much perturbations of the azimuthal component 
of the magnetic fields (cf.. Figs. 5-6). Note that in the linear 
limit, Vq is governed by the wave equation of Eq. m which 
was derived by Hollweg (1981). This equation differs from 
Eq. (|32]) . and this difference affects the evolution scenario 
of velocity perturbations. 
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Fig. 6. Spatial profiles of Bz{x,y,z = 0) associated with torsional Alfven waves, at t = 50 s (top-left), t = 100 s (top-right), 
t = 150 s (bottom-left), and t = 200 s (bottom-right). A part of the simulation region is displayed only. 


Torsional Alfven waves can be traced in the vertical 
profiles of Vz{x^y,z = 0) = Vo{x^y^z = 0) (see Fig. [TQ|) . 
At t = 200 s (the bottom-right panel), the upwardly prop¬ 
agating Alfven waves resulting from the initial pulse are 
well seen, with their amplitudes of |14| ~ 16 km s“^. The 
downwardly propagating waves subside as it is discernible 
at t = 100 s (the top-right panel) and at t = 150 s (the 
bottom-left panel). At t = 100 s, the upwardly propagating 
waves reach the altitude of more than ^ = 1.6 Mm as shown 
in the top-right panel of the figure. 


As a result of nonlinearity, torsional Alfven waves drive 
a vertical flow through a ponderomotive force (Hollweg et 
al. 1982; Nakariakov et al. 1997, 1998). We refer to Eq. (6) 
in Hollweg et al. (1982) which illustrates all the terms which 
contribute to the component of velocity that is parallel to 
a magnetic field line, . Denoting by s distance measured 
along this line and by r{s) a distance of any point on this 
line from the ^-axis we restate Eq. (6) in Hollweg et al. 
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Fig. 8. Maximum of Vo (left-panel) and Be (right-panel) vs height. The numerical (analytical) data points are represented by 
asterisks (solid line). 
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Fig. 7. The top panel: ratio of magnetic to kinetic energies of 
Alfven waves; the bottom panel: ratio of slow magnetoacoustic 
waves kinetic energy to Alfven waves energy (asterisks) and ratio 
of fast magnetoacoustic waves kinetic energy to total energy 
of Alfven waves (diamonds). These energies are evaluated by 
integration in space over the entire computational domain. 


(1982) as 

^ f , A (Al) = _J_^ + ^ + 

dt \ B|| ) ds \ B|| ) B|| 5s B|| 


1 


qV^ 


1 2 ^ d\nr 


ds\2ii) 


(34) 


where g\\ and are the components of, respectively, 
the gravitational acceleration and perturbed magnetic field 
along the magnetic field line, and Vq and Bq are azimuthal 
components of velocity and magnetic field, respectively. 
Following Hollweg et al. (1982), we state that the terms 
in the square bracket result from the s-component of the 
Lorentz force: the first term denotes the 5 -component of 
the centrifugal force that corresponds to twisting motions; 
the second term represents the magnetic tension; the third 
term is the magnetic pressure. 

The vertical flow can be seen in the spatial profiles of 
Vy{x, y,z = 0) (Fig. [TTl middle). At t = 400 s, slow magne¬ 
toacoustic waves in a low plasma /3 region of the magnetic 
fiux tube of axial symmetry as displayed in the right-panel 
of Fig. [3] are weakly coupled to fast magnetoacoustic waves, 
which are described by essentially only Vy component of 
the velocity. Our simulations demonstrate that at t = 175 s 
slow magnetoacoustic waves already penetrated the solar 
corona, reaching the level higher than ^ = 4.5 Mm. 

It is noteworthy that non-linear torsional Alfven waves 
drive perturbations in a mass density, which are associated 
with fast (Fig. [TTl f^op) predominantly slow (Fig. [TTl 
middle) magnetoacoustic waves. The vertical spatial profile 
of log(p(x, z = 0)1 go) is shown at t = 225 s in the bottom 
panel of Fig. [TTl Here go = 10 “^^ kg m“^ is the hydrostatic 
mass density at the level ^ = 10 Mm. At this time the dense 
plasma jets reach the altitude of ^ ~ 3.5 Mm. The plasma is 
lifted up due to the magnetic pressure gradient originating 
from the Alfven waves and acts against the gravity to push 
the chromospheric material up to the inner solar corona. 
Therefore, such predominant vertical pressure can create 
the cool plasma jet flows typically observed in the solar 
corona. 

Figure [ 12 ] displays the time-signatures of T4 (the solid 
line) and Vy (the dashed line) collected at the detection 
point, {x = 0 . 2 ,^ = 3 . 5 , 2 ) = 0 . 2 ) Mm. We infer that slow 
magnetoacoustic waves, represented by Vy, reach the detec¬ 
tion point at t ^ 130 s, while Alfven waves arrive at this 
point about 30 s earlier. 

Noting that up to y ^ 0.7 Mm Alfven speed is lower 
than the sound speed (Fig. [3] left), the travel times of the 
torsional Alfven, tAt(^), and slow magnetoacoustic waves. 
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Fig. 9. Streamlines at t = 100 s (top-left), t = 150 s (top-right), t = 200 s (bottom-left), and t = 225 s (bottom-right). Red, green 
and blue arrows design the directions of the x-, y-, and z-axis, respectively. The size of the box shown is (—0.75, 0.75) x (0,6) X 
(—0.75, 0.75) Mm. The color map corresponds to the magnitude of a total velocity. 


tst{y), are given by 



dy 

CA(r, y) 
dy 

CT{r,y) 


(35) 

(36) 


where the integration is performed along the corresponding 
magnetic field line, and ct = CsCa/\/c^ + c\ is the tube 
speed. Using Eqs. dsn and (|36]), we find tAt(^ = 1-9) ~ 
157.7 s and tst{y = 1.9) ^ 173.3 s. These values are larger 


than the arrival times of Alfven and slow magnetoacoustic 
waves. This dispatch results from nonlinearity. 

The amplitude of the Alfven waves, seen in V^, pene¬ 
trating into the upper region of the atmosphere drops from 
about Vo = 9 km s“^att = 100s just under the transition 
region, through Vq = 8.4 km s“^ just above the transition 
region. We can evaluate the transmission coefficient. 



(37) 
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Fig. 10. Spatial profiles of Vz{x,y,z = 0) associated with torsional Alfven waves at t = 50 s (the top-left panel), t = 100 s (the 
top-right panel), t = 150 s (the bottom-left panel) and t = 200 s (the bottom-right panel). 


where Ai (At) is the amplitude of the incident (transmitted) 
waves. Substituting Ai = 9 km s“^ and At = 8.4 km s“^ 
into the above formula we get Ct ~ 0.95, which means that 
about than 5% of the waves amplitude became reflected in 
the transition region and about 95% was transmitted into 
upper atmospheric layers. These are rough estimates due 
to difficulties with precise estimation of wave amplitudes 
which highly vary in time at the transition region. 


4. Discussion and Conclusions 

In this paper, we present an analytical model of a mag¬ 
netic flux tube of axial symmetry embedded in the solar 
atmosphere whose photosphere, chromosphere and transi¬ 
tion region are described by the model of Avrett & Loeser 
(2008). Our flux tube model can easily be adopted to any 
axisymmetric magnetic structure by specifying a different 
magnetic flux. The model was used to perform numerical 
simulations of the propagation of torsional Alfven waves 
and their coupling to magnetoacoustic waves by using the 
FLASH code. Torsional Alfven waves are launched at the 
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Fig. 12. Time-signatures of velocities Vz (solid line) and Vy 
(dashed line), collected just below the transition region at the 
detection point {x = 0.1, y = 1.9, z = 0.1) Mm. 


top of the solar photosphere by an initial localized pulse in 
the azimuthal component of velocity. This pulse mimics the 
convectively excited photospheric vortices or any other kind 
of plasma rotatory motions recently observed in the solar 
atmosphere (see Section 1 for references and discussion). 

Our results show a complex behavior of torsional Alfven 
waves and the driven magnetoacoustic waves, and their re¬ 
sponses in the solar chromosphere, transition region, and in¬ 
ner corona. The time and spatial evolution of Alfven waves 
shown by our numerical results is akin to the dynamics of 
fast solar swirling motions observed in the form of various 
jets at diverse spatio-temporal scales. Therefore, our model 
presented here has a potential to describe the physics of 
various observed and recurring swirling motions (jets) in 
the solar atmosphere in which the driver excites torsional 
Alfven waves as well as the associated plasma perturbations 
(e.g., Wedemeyer-Bohm et al. 2012; Murawski et al. 2014, 
and references cited therein). 

Since in our model the solar magnetic flux tube of 
axial symmetry is strongly magnetized in the solar pho¬ 
tosphere, it can be either associated with strong mag¬ 
netic bright points (MBPs) (Jess et al. 2009) or with the 
concentration of the magnetic field between intergranular 
lanes (Wedemeyer-Bohm et al. 2012). With the use of the 
Swedish Solar Telescope Jess et al. (2009) observed the pe¬ 
riodic variations above such MBPs in the detected (at dif¬ 
ferent altitudes) Doppler widths of the line. The authors 
claimed that their findings reveal the presence of the prop¬ 
agating torsional Alfven waves with their wave periods of 
2 — 12 min up to the top of the solar chromosphere in the 
expanding wave-guide rooted in the strongly magnetized 
MBP. Clearly, our developed model can be used to excite 
torsional Alfven waves and comparing the theoretically pre¬ 
dicted range of wave frequencies to that established obser- 
vationally. Indeed, this was one of the main goals achieved 
in our paper and we found the theoretical wave periods 
to be smaller than 100 s (Fig. [12)), which are close to the 
observational findings of Jess at al. (2009). 

Our numerical results also show that the excited per¬ 
turbations in azimuthal component of the magnetic field 
(Fig. [6]) are seen below the solar transition region. However, 
in the inner corona these perturbations are of a significantly 


Fig. 11. Spatial profiles of V^(x^y^z = 0) at t = 100 s (top), 
Vy{x,y,z = 0) at t = 175 s (middle) and \og{g{x,y,z = 0)/^o) at 
t = 225 s (bottom), where go = 10“^^ kg m“^ is the hydrostatic 
mass density at the level y = 10 Mm. 
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lower amplitude. On the other hand, the perturbations in 
azimuthal component of velocity are stronger in the solar 
corona. 

Finally, let us point out that our analytical model can 
easily be adopted to derive the equilibrium conditions for 
any axisymmetric 3D magnetic structure. Our devised nu¬ 
merical model of the solar magnetic flux tube of axial sym¬ 
metry can be applicable to the variety of flux tubes with 
different photospheric fields, and different expansions, as 
well as with different strengths of the drivers. The model 
can also be used to investigate the nonlinear evolution of 
coupled torsional Alfven and magnetoacoustic waves, and 
to determine and understand drivers of coronal jets. 
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